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We study theoretically the surface of a topological insulator with a sliding magnetic superlattice 
coated above. By analyzing time-dependent Dirac equations, the dynamics of the zero mode is 
investigated. When the superlattice's sliding velocity is smaller (larger) than the Fermi velocity 
of topological insulator, the zero mode is perfectly (imperfectly) pumped. We also propose the 
application of this setup, which rectifies currents or generates pulse currents. It would provide a 
prototype of electronic devices based on topological insulator. 



Introduction: Topological insulator [lMJl is a new state 
of matter with insulating bulk and metallic surface. The 
gapless surface is protected by the topology, and persist- 
ing even in the presence of disorder as long as time rever- 
sal symmetry is respected. Both theory and experiment 
have established that this surface is not a conventional 
metal, but a helical liquid. As a result, a dissipationless 
spin current can be realized on the surface. 

When time reversal symmetry is broken on the surface, 
even richer properties of topological insulator emerges, 
such as monopoles[5] and chiral Majorana fermions[6j. 
When a magnetic domain wall is attached on the sur- 
face, a chiral mode is generated inside the wall, which 
shares the same physics as Jackiw and Rebbi's zero 
modeQ- In this sense, topological insulator bridges be- 
tween condensed matter physics and high energy physics. 
Many unrealized phenomenon in high energy physics are 
promisingly addressed experimentally in topological in- 
sulator. 

However, potential applications of topological insula- 
tor are seldom mentioned. Let us consider a magnetic 
domain-wall lattice created on its surface. Then, do- 
main walls generate series of chiral zero modes, which 
are treated as quantum wires. By the motion of these 
domains, we can tune the electronic signals. It would 
provide a prototype of application in electronics. 

In this paper, we find a dissipationless transverse cur- 
rent induced by a magnetic superlattice attached onto 
a topological insulator by solving time-dependent Dirac 
equations. The behavior changes drastically whether the 
sliding velocity of the superlattice is larger or smaller 
than the velocity of the Dirac particle. The transverse 
current is perfectly pumped when we apply a slowly mov- 
ing magnetic lattice (adiabatic quantum pumping) and 
imperfectly pumped when we apply a fast moving lat- 
tice (non-adiabatic quantum pumping). We also sug- 
gest an application of this set up to a current modulator, 
which converts between direct, alternating, and repulse 
currents. 

Model: Our system consists of a three-dimensional 
topological insulator with a ferromagnet attached onto 
it [FigH]. Electrons obey 2D Dirac equations on the 




FIG. 1: Illustration of an experimental setup of current mod- 
ulator on the surface of a topological intulator. 



surface of a topological insulator, as has been demon- 
strated by the spin- and angle-resolved photoemission 
spectroscopy [8J. In the presence of a ferromagnetic layer 
on top of a topological insulator, the Hamiltonian is given 
by 

H = hv F (k x cr x + kyCF y ) + /i B M (x, t) ■ a, (1) 

where is the Fermi velocity of Dirac fermions, which 
is 3 x 10 5 m/s for Bi2Se3[9|, oi are the Pauli matrices, 
Hb is the Bohr magneton, and M (x, t) represents the 
strength of the exchange coupling. It has been argued[10j 
that the in-plane magnetic field contributes only tiny ef- 
fects to the transport property of a topological insulator. 
Namely, it is a good approximation that M (x, t) has only 
the z-component, M (x, t) = (0, 0, M (x, t)). Let us as- 
sume that the system is homogeneous along the ?/-axis. 
Then, the central point is the x and t dependence of this 
exchange coupling, which is realized by domain wall mo- 
tion along the x-axis with velocity v in the ferromagnet 
layer. It is reasonable to set 



M(x,t) 



. [(x - vt) /A] 



(2) 



to describe a domain- wall lattice, where m is a numerical 
constant characterizing the strength of a single domain. 
We define lj = v/X. 

Static solution: The sinusoidal behavior of M{x,t) 
leads to a series of zero-field points periodically. As 
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FIG. 2: We illustrate the magenetization due to the mag- 
netic superlattice M (x) (dashed curve) and the static solu- 
tions \ip (x)\ 2 associated with it (solid curve), (a) for the + 
sign and (b) for the — sign in (|6)). The horizontal line repre- 
sents the x axis. We have set /jbttiX/ '(Hvf) = 5. 



particle-hole symmetry is respected, according to Jackiw 
and Rebbi's work, zero modes appear around these 
points. In order to derive the explicit expression of the 
zero modes here, let's start from the time-dependent 
Schrodinger equation 



(3) 



where we have set the wave function as \I/ = (^,^) . 
We may set k y =constant in (pQ) due to the translational 
invar iance along the y- axis. 

We first study the static case where uo = 0. Particle- 
hole symmetry guarantees the existence of zero-energy 
solutions with the relation x/j = ^ = ±£0j,. When k y = 
0, the equation of motion is transformed into 



We may solve it asQ 



± 



ip± (x) = exp 
which yields 

ijj± (x) = exp 



f 

JO 



M(x')dx' 



^/i^mA x 
hvp X 



(5) 



(6) 



up to a normalization constant. We illustrate the mag- 
netization and (x)\ 2 for these two cases in Figj2j 

We may also present solutions for the non-zero k case. 
The two solutions have opposite chirality, so their group 
velocity along ^-direction is opposite as well. The wave 
function is a linear combination of ip± (x), 



(7) 



with ©. 

Moving Domain Wall: We start with the investigation 
of the domain wall motion described by ([2]). The sinu- 
soidal potential has a double-periodicity, 

M (x + L, t) = M (x, t), M(x,t + T)=M (x, t) , (8) 

where L = 2ttX and T = 2ttX/v. To such a system 
both the Bloch theorem and the Floquet theorem are 
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FIG. 3: We show numerical solutions of the time-dependent 
Dirac equation with the sinusoidal potential for (a,a') uj/uj c = 
0.25, (b,b') uj/uj c = 3, (c,c') uj/uj c = 10. The vertical axis is 
\ip (x, t)| 2 . The horizontal axes are x and t. (a), (b) and (c) are 
bird's eye view and (a'), (b') and (c') are their contour plot. 
We also present analytic solutions in (a") and (c"), which are 
valid cj/cj c <C 1 and u/u c ^> 1, respectively. The figures are 
plotted in the region 0<x<L,0<£<T for the adiabatic 
case, and in the region < x < L, < t < T c for the non- 
adiabatic case, where T c = (uj/uj c )T. We have started from 
the static solution with fismX/ '(Hvf) = 1. 



(4) applicable: The wave function is of the form ijj a (x, t) 



e l£oit+lkoiX u ot (x,i) with u a (x + L,i) = u a (x,t) and 
u a (x,t) = u a (x,t + T). Then it is enough to analyze 
the torus region < x < L,0 < t < T. 

First we investigate the case where the velocity of the 
moving lattice is low, v < i>f, or uo < uo c with uo c = vp/X. 
We have carried out a numerical analysis of the time- 
dependent Dirac equation, starting from a static solution, 
whose results we give in Figj3^a,a'), where we have set 
uj/uj c — 0.25. We find that the zero mode moves together 
with the magnetic lattice, and it is adiabatically pumped. 

It is possible to derive the analytic solution based on 
the adiabatic approximation, 



ip± (x, t) = exp 



/iB^A (x — vt) 
±— cos 



X 



(9) 



which is constructed by making the Galilei boost of the 
static solution ([5]): See Fig|3l(a"). This solution is valid 
when the velocity is very low. There are fluctuations in 
\ip {x, t) | in the case of the numerical calculation, which is 
absent in the analytic one. This is because the adiabatic 
solution ([9]) is not an exact solution of (|3j). In any case, if 
we introduce the function P n (t) representing the position 
of \iJj (x, t)|'s n-th peak at time t in unit of A, we have 

Pn (T) ~ P n (0) = 1 (10) 

for uo < uj c . In that sense, the zero mode is adiabatically 
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pumped along x-direction as long as the velocity of do- 
main wall motion is smaller than the Fermi velocity. As 
the Fermi velocity for topological insulator is generally 
very large, this requirement is usually satisfied. Due to 
the periodicity along x-direction, P n (t) actually doesn't 
depend on n , and the subscript is neglected in the fol- 
lowing. 

It's also interesting to investigate the case where the 
velocity of the moving lattice is fast, v > vp, or uj > 
uj c . We have carried out a numerical analysis, whose 
results we give in Figj3fb,b') and (c,c'), where we have 
set uj/uj c = 3 and 10. Fig[3](b,b') are plotted in the region 
< t < T, while Figj3](c,c') are plotted in the region 
< t < T c with T c = (uj/uj c )T. 

In the large uj limit, we give also an analytic solution. 
Since the magnetization oscillates very quickly, we can 
approximate the system by the one-period time-averaged 
equation of motion(l3|. 

ikd t ^ = H^, (11) 

where 

H=^j^ H(x,t)dt. (12) 

The equation of motion is rewritten as 

ihdt^ = -ihv F (d x a x + dy(jy) (13) 

or 

d 2 t ^ = vl{d 2 x + dl)^. (14) 

It does not depend on uj. The solution is given by 
[ip (x — v-pt) + ip (x + vpt)] /2 with ([6]), which we display 
in Fig|3](c"). The agreement is good between the ana- 
lytic and numerical solutions when the velocity is large 
(uj/uj c > 10). 

The pumping rate is independent of a;. It implies that 
one transverse current is pumped in the period which is 
not T but T c , 

P(T C )-P(0) = 1. (15) 

The primitive torus < t < T has been extended to 
< t < T c . As a result, the pumping transverse current 
in one period T is anti-proportional to the frequency, 

P(T)-P(p)=u c /u. (16) 

The reason of this imperfect pumping is that the trans- 
verse current can not follow the magnetic lattice because 
it moves too fast. 

A comment is in order. The adiabatic pumping case is 
experimentally realsitic, since it is hard to move domain 
walls faster than v > v-p experimentally at this stage. 
Nevertheless the non-adiabatic case is theoretically very 
interesting, as we have shown. 



(a) Adiabatic pumping (b) Non-adiabatic pumping 
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FIG. 4: Schematic diagram of the adiabatic (a) and non- 
adiabatic (b) pumping. One transverse current is pumped in 
the time T in the adiabatic pumping, while one transverse 
current is pumped in the time T c in the non-adiabatic pump- 
ing. Numerical results for the non-adiabatic pumping for the 
period (c) [0,T], and (d) [0,T C ]. 




FIG. 5: u dependence of P (T) — P (0). The horizontal axis 
is uj and the vertical axis is P (T) — P (0) . 

In order to analyze the transition between the adi- 
abatic and non-adiabatic pumpings, we calculate the 
pumping transverse current numerically. The pumping 
rate is estimated from the velocity of the position where 
|^ (x,t)\ takes the maximum value at t = T. We show 
the numerical solution in FigJS] It is seen that the pump- 
ing transverse current is well described by 

p v- p «»={^ ! z :>>::■ <"> 

The analytic and numerical solutions show a good 
agreement. The change between the adiabatic and non- 
adiabatic transverse current pumping takes place rather 
suddenly. It may be a kind of a non-equiribrium phase 
transition. 

Experimental setup: The realization of magnetic su- 
perlattice is a highly nontrivial question. Here we 
present an experimental proposal designed in the adia- 



batic pumping regime, as illustrated in FigJTJ On top 
of a topological insulator, a ferromagnetic thin layer is 
coated. Above the ferromagnetic layer, an array of peri- 
odically distributed superconductors is deposited. Usu- 
ally there are many domains in the ferromagnet, which 
are naturally generated and hard to control. However, 
with the help of this array of superconductors, it be- 
comes possible. At the beginning, an upward external 
magnetic field magnetizes the whole ferromagnetic layer 
in the same direction. Then, we lower the temperature 
below the critical temperature of superconductor, and re- 
verse the direction of external magnetic field. Due to the 
Meissner effect, magnetic field is screened just below the 
superconductor, and unscreened elsewhere. As a result, 
magnetization in the areas without superconductor above 
is reversed, and the magnetization below the supercon- 
ductor may remain unchanged if the external magnetic 
field is properly controlled. We would obtain a magnetic 
superlattice on top of a topological insulator in this way. 

It is well known [13] that a current can drive the mag- 
netic motion of domain walls. By introducing a current in 
the ferromagnetic layer along the x-direction, we can ap- 
proximately realize the moving magnetic configuration in 
([2]), and conducting channels provided by the zero modes 
are generated. Once the device is bridged by a conduct- 
ing channel, a significant current is detected. 

It's an interesting property that the two zero modes 
in one period of the magnetic superlattice have opposite 
chiralities. These two conducting channels have tendency 
to transport opposite currents. We propose three type 
devices. 

(a) When we do not apply the current into the ferro- 
magnet, the magnetic superlattice is static. We apply an 
alternating voltage between two leads parallel to the do- 
main wall. The resulting current is rectified into a direct 
current, 



input 



output 



Jout = max [J in ,0] , 



(18) 



because the zero-energy conducting channel is chiral and 
hence only the forward bias current goes through. This 
system rectifies alternating currents to direct currents, or 
acts as a diode. 

(b) When we apply the current into the ferromagnet, 
the domain wall moves[13j. We apply the direct voltage 
between two leads. The resulting current is modulated 
to be a pulse-current, 



j+(t) 



:Re 

fliky 

2m 



h 
2mi 



exp 



llbttiX (x — vt) 
2 — cos 



X 



(19) 



because the zero energy conducting channel passes 
through the contact periodically according to the motion 
of domain walls. This system acts as a pulse generator. 
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FIG. 6: We show the time dependence of transverse current, 
(a) current rectifier (b) pulse generator (c) alternating pulse 
generator. We have set jiBrnX/ (Kv?) = 10. 



(c) When we apply an alternating voltage to the same 
setup with (b), the resulting signal between two devices 
is an alternating pulse. This current is given by 



J (t) = Re 

hk 



2mi y 



- sinh 



_ /ieraA (x — vt) 
2— cos 



X 



(20) 



and is shown in Figj6] This system acts as an alternating 
pulse generator. 

It will be possible to measure these currents by at- 
taching leads to the topological insulator parallel to the 
domain wall as in FigQ] 

Conclusion: In this paper, we have studied the dy- 
namics of the zero mode in the presence of a magnetic 
superlattice on top of a topological insulator, which is de- 
scribed by a time-dependent Dirac equation. The domain 
wall motion of the superlattice is studied. We have pro- 
posed a prototype of electronic device based on the this 
theoretical studies, where the input current is found to 
be significantly modulated by the magnetic superlattice 
[FigH] . We hope this work can promote the application 
of topological insulator in future. 
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